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Abstract 

c/3 ' Given a planar undirected n-vertex graph G with non-negative 



edge weights, we show how to compute, for given vertices s and t 
in G, a min st-cut in G in 0(n log log n) time and 0{n) space. The 
previous best time bound was O(rilogn). 
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^ ■ 1 Introduction 

en ■ 

Q . Given a graph G = (V, E) with non-negative edge weights and given vertices 

O ! s,t & V, an st-cut of G is a partition of V into two subsets S and T = V\S 

such that s ^ S and t E T. The weight of the st-cut is the sum of weights 
of edges starting in S and ending in T. A min st-cut of G is an st-cut of G 
^ ■ having minimum weight. 

Computing a min st-cut is a classical algorithmic problem with several 
applications in areas such as chip design, communication networks, trans- 
portation, and cluster analysis. The problem is intimately related to another 
well-studied problem, that of computing a max st-fiow. The classical max 
flow min cut theorem implies that the weight of a min st-cut is the value of 
a max st-flow. 

For general graphs, the fastest known max st-flow algorithm runs in time 
0{'mn\og{m? /n)), where m is the number of edges and n is the number of 
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vertices [7] . Another implication of the max flow min cut theorem is that a 
min st-cut can be obtained from a max st-flow in linear time. Hence, a min 
st-cut can also be computed in Olmnloglrn^/n)) time. 

For planar undirected graphs, Reif |Tl] showed how to solve the min st- 
cut problem in 0(?T,log^n) time. This was later improved to O(nlogn) by 
Frederickson [B]. For directed planar graphs, Borradaile and Klein [2] showed 
that the max st-fiow problem can be solved in 0(?7,log?7,) time and this gives 
an O(nlogn) time min st-cut algorithm for such graphs. 

In this paper, we give a min st-cut algorithm for planar undirected graph 
with O(nloglogn) running time and 0{n) space requirement, thereby im- 
proving the time bound by Frederickson [B]. In order to achieve this, we do 
not depart from Reif's approach [TT]. Instead we speed it up using a two- 
phase approach. The first phase runs a "coarse" version of Reif's algorithm 
which only determines a subset of the min st-cut candidates found by the 
original algorithm. We obtain a running time of O(nloglogn) for this phase 
using the fast Dijkstra variant of Fakcharoenphol and Rao [1]. In the second 
phase, the remaining min st-cut candidates are found exactly as in the algo- 
rithm by Reif but since the first phase partitions the problem into simpler 
subproblems, we can show that the second phase also runs in O(nloglogn) 
time. 

The organization of the paper is as follows. In Section |2l we give some def- 
initions and introduce some of the tools that we need. We briefly go through 
the ideas of Reif's algorithm in Section [3] before presenting our algorithm 
in Section |H One step of our algorithm constructs a certain division of the 
graph and we present the details of this step in Section |5l Finally, we make 
some concluding remarks and suggestions for future research in Section |6l 

2 Preliminaries 

For a graph G = (V, E), define a piece P = {Vp, Ep) of G to be the subgraph 
of G defined by a subset Ep of E. In G, the vertices of Vp incident to vertices 
in V\Vp are the boundary vertices of P. Vertices of Vp that are not boundary 
vertices of P are interior vertices of P. If G is edge-weighted, we define the 
dense distance graph of P to be the complete graph on the set of boundary 
vertices of P where each edge (m, v) has weight equal to the shortest path 
distance (w.r.t. the edge weights) in P between u and v. 

Let G = (y, E) be an n- vertex planar graph with a non-negative weight 



function w : ^ — ?> M defined on its vertices. For a subset A of V, define 
w{A) = S^e^w(t')- We assume that w{V) = 1. The separator theorem of 
Lipton and Tarjan states that in 0{n) time, V can be partitioned into three 
subsets A, B, and C such that 

• no edge joins a vertex in A with a vertex in B, 

• I < w{A) < I and I < w{B) < |, and 

Using this theorem, Frederickson [6] showed how to obtain, for any pa- 
rameter r G (0,n), an r-division of G, which is a division of (the edges of) 
G into 0{n/r) pieces each containing 0{r) vertices and 0{y/r) boundary 
vertices. He gave an 0(nlogr + {n/y/r) logn) time algorithm to find such a 
division. 

A stronger version of the separator theorem is the cycle separator theorem 
of Miller [10] which states that if G is a plane graph then C can be chosen 
such that there exists a Jordan curve that only intersects G in vertices of C. 
Miller showed that such a separator can be found in linear time. 

We will show that by applying the cycle separator theorem as well as ideas 
of Fakcharoenphol and Rao [4], we can obtain the r-division of Frederickson 
but with some additional properties. More precisely, define the holes of a 
piece to be the internal faces containing boundary vertices. We prove the 
following result in Section |5l 

Theorem 1. For a plane n-vertex graph, an r-division in which each piece 
has 0(1) holes can be found in 0{nlogr + (n/y/r) logn) time. 

In the following, when we talk about an r-division, we shall assume that 
it has the form in Theorem [H 

We shall identify an st-cut with the set of edges from the s-side to the 
t-side of the cut. 



3 Reif s Algorithm 

Reif's algorithm [TT] makes use of the following duality between cuts in the 
primal graph and cycles in the dual graph: a min st-cut in a plane undirected 
graph G corresponds to a minimum weight simple cycle separating faces s 




Figure 1: (a): In cut-open graph G*^, Reif s algorithm computes a shortest 
path Q from the midpoint on P to the midpoint on P' and recurses on 
the two subgraphs generated, (b): The coarse version of Reif's algorithm 
only computes shortest paths between boundary vertices on the cut-path. A 
refined version is then applied to find the remaining shortest paths. Only 
shortest paths from the coarse version are shown. Dashed line segments show 
the boundaries of pieces in the r-division. 



and t in the dual G* of G; here, a simple cycle is said to separate two faces 
if one face is in the interior and the other is in the exterior of the cycle. 

In the first step of Reif's algorithm, a shortest path P = pi — )■ p2 ^ 
■■■—)■ p\p\ from an arbitrary vertex pi on face s to an arbitrary vertex p|p| on 
face t in G* is computed. Then G* is "cut open" along P as follows. Remove 
the set Er of edges emanating right of P in the direction from s to t. Insert 
a copy P' = p']^ — 7- P2 —!■■■■ —)■ P|p| of P and for each edge {pi, u) E E^, add 
edge {p'i.u). We let G^^ be the resulting graph, see Figure [U^a). 

Next, Reif's algorithm computes a shortest path Q in G*^ from the mid- 
point p\\p\/2\ of P to the midpoint Pnp|/2i of P' ■ This splits G*^ into two 
subgraphs and splits P and P' into two halves, one for each side of Q. In 
each of the two subgraphs, degree two- vertices are removed by merging their 
incident edges. The algorithm then recurses on the two subgraphs and the 
two subpaths. 

Let Qi be the shortest path found by the algorithm and let pi and p[ be 
the first and last vertex of Qi, respectively. Then the cycle in G* obtained 
from Qi by identifying pi with p[ is a minimum-weight st-separating cycle 
in G* . By the min cycle/min cut duality, this cycle defines a min st-cut in 
primal graph G. 



With Dijkstra's shortest path algorithm, Reif s algorithm runs in 0(n log^ n) 
time. This can be improved to O(nlogn) time with Frederickson's algo- 
rithm ^ or by speeding up Reif s algorithm using the linear time shortest 
path algorithm of Henzinger et al. [8]. In the next section, we will further 
improve Reif's algorithm to get O(nloglogn) running time. 

4 The O(nloglogn) Time Algorithm 

In this section, we present our algorithm and give the claimed O(nloglogn) 
time bound. To ease the presentation, we leave out some details of the 
algorithm and return to them in Section 14. 3[ 
We start with the following simple lemma. 

Lemma 1. Let s and t be faces in a planar n-vertex graph G and let P 
he a given shortest path between a vertex on s and a vertex on t. In an 
application of Reif s algorithm to find a minimum weight st-separating cycle 
in G, consider a subproblem defined by a subgraph H and a subpath of P 
of length 0{\og^n) for a constant c. Then this subproblem can be solved in 
0(|i/| loglogn) time. 

Proof. Recursion depth for Reif's algorithm in H is only 0(log(log'^?7,)) = 
O (log log n) so the running time for the subproblem is 0{\H\ log log n) using 
the linear time shortest path algorithm in [8]. D 

We essentially run Reif's algorithm but speed part of it up with the 
Dijkstra variant of Fakcharoenphol and Rao |1]. In the following, let G 
denote the dual of a plane embedding of the input graph. We need to find a 
minimum weight st-separating cycle in G, where s and t are faces. 

Let P be a shortest path in G from an arbitrary vertex on s to an arbitrary 
vertex on t. We can find this path in linear time using the algorithm in [8] . We 
first run a "coarse" version of Reif. This will identify in O(nloglogn) time a 
subset of all the st-separating cycles found by the original algorithm. More 
precisely, the set of cycles found will split G into subgraphs each of which 
contains a subpath of P of length 0(log'^ n) for some constant c. We then run 
the "refined" Reif algorithm by applying Lemma [1] to each subgraph and its 
associated subpath of P. This will find the minimum weight st-separating 
cycle in G. By ensuring that the total size of the subgraphs is 0{n)., the entire 
algorithm runs in O(nloglogn) time. Figure [It^b) illustrates the output of 
the first phase of our algorithm. 



4.1 First phase 

We will now describe the first phase of our algorithm which is the coarse 
version of Reif 's algorithm. 

r-division First, we apply Theorem [1] to obtain an r-division of G for 
r = log n. This takes 0{n\ogr + [n/ y/r) logn) = 0(?T,loglogn) time. 

Cutting pieces open We make an incision in G along shortest path P as 
in Reif 's algorithm. This induces incisions in those pieces containing parts 
of P and we update the pieces accordingly. If a boundary vertex of a piece 
belongs to P before the incision, we regard both of its two copies after the 
incision as boundary vertices of that piece. Note that there will still be only 
0{^/r) boundary vertices in each piece and these boundary vertices will still 
be on a constant number of faces after the incision. Hence, the resulting set 
of pieces forms an r-division in the cut-open graph. 

Dense distance graphs Next, we compute dense distance graphs of the 
pieces in the r-division. To do this, we shall apply Klein's multiple-source 
shortest paths algorithm [U] . For an /i- vertex plane graph H and a fixed face 
/ of iJ, this algorithm builds a data structure in 0{hlogh) time and space 
such that shortest path queries between vertex pairs {u, v), where either u or 
V is on /, can be answered in 0{\ogh) time per query. 

For each piece, we apply Klein's algorithm to set up a data structure 
for the external face and query this data structure for shortest path dis- 
tances in the piece from boundary vertices on this face to all other bound- 
ary vertices. Since there are 0{r) pairs of boundary vertices, we obtain 
all these distances in 0(r log r) time. We similarly set up data structures 
for each hole and get the distances between the remaining pairs of bound- 
ary vertices. Since the piece has a constant number of holes, total time 
to construct its dense distance graph is O(rlogr). Over all pieces, this is 
0((r2/r)r log r) = 0{n\ogr) = O (n log log n) time. We represent the edge 
weights of each dense distance graph in a distance matrix with 0{^/r) rows 
and columns. 

Fast Dijkstra Fakcharoenphol and Rao [1] gave an efficient Dijkstra vari- 
ant for planar graphs. More precisely, they showed the following. Given a 



collection of pieces each having a constant number of holes and given their 
dense distance graphs, any shortest path tree in the union of these dense 
distance graphs can be computed in 0{hlog n) time, where h is the total 
number of vertices in these graphs. In general, this time bound is sublinear 
in the number of edges. 

Since we have computed dense distance graphs for the pieces in our r- 
division and since the total number of boundary vertices of these pieces is 
0{n/^/r), it follows that a shortest path between any two boundary vertices 
in the r-division can be computed in 0{{n/y/r)log n) = 0{n/ logn) time. 
Note that this shortest path consists of edges from the dense distance graphs 
so it is an implicit representation of a shortest path in the underlying cut- 
open graph. 

Coarse Reif Let pi, . . . , p\p\ be the ordered sequence of vertices of P, start- 
ing with the vertex on face s. Let Pi^, ■ ■ ■ ,Pif. be the (possibly empty) sub- 
sequence of vertices that are boundary vertices in pieces of the r-division. 
Note that these vertices partition P into subpaths each of which is contained 
in a piece. 

The coarse version of Reif's algorithm is the normal algorithm of Reif 
restricted to subsequence Pi^, . . . ,Pi^. and using the 0(n/logn) time shortest 
path algorithm for each vertex in the subsequence. Since recursion depth is 
O(logn), total time for the first phase of our algorithm is 0{n) in addition 
to the 0(n log logn) time to find the r-division and to set up the dense 
distance graphs. The st-separating cycles found in this phase partition G 
into subgraphs each containing a subpath of P fully contained in a piece of 
the r-division. Hence, the length of each such subpath is bounded by the 
size 0{r) = 0{\og^ n) of a piece. 

Subgraphs for recursive calls When applying the coarse version of Reif's 
algorithm, we need to find the subgraphs for recursive calls. Consider a sub- 
graph H in some recursive call. We associate with H the boundary vertices 
belonging to H and the cyclic orderings of these vertices on holes and ex- 
ternal faces of pieces. When running the fast Dijkstra variant for H, only 
distances between boundary vertices associated with H are considered in the 
dense distance graphs. Hence, these graphs need not be updated in recursive 
calls. The time to find a shortest path in H will be 0{hlog n), where h is 
the number of boundary vertices in H. 



4.2 Second phase 

In order to run the second phase of our algorithm, we need to convert the 
shortest paths consisting of edges from dense distance graphs to the underly- 
ing shortest paths in G and we need to find the subgraphs of G bounded by 
these paths. In Section U3] we show how to do this in 0{n) time such that 
the total size of the subgraphs is 0{n). Applying Lemma [T] with constant 
c = 6 to each subgraph, we get O(nloglogn) time for the second phase of 
our algorithm. Hence, the entire algorithm has 0(n log log n) running time. 

4.3 Overlapping subgraphs 

We need to ensure that the total size of the subgraphs generated by the 
coarse version of Reif 's algorithm is not too large, i.e., we need the total size 
to be linear. The problem is that subgraphs overlap so a vertex can belong to 
several subgraphs. The original algorithm of Reif ensures linear total size by 
deleting, in every subgraph generated, each degree two vertex by replacing 
the two edges ei and 62 incident to it by one whose weight is the sum of the 
weights of ei and 62- 

First phase In the first phase of our algorithm, we do something similar: 
if a boundary vertex p is incident (in the union of dense distance graphs) 
to only two other boundary vertices gi and g2 for the current subgraph, p is 
removed and a single super edge is added between gi and ^2 whose weight is 
equal to the sum of weights of edges (gi, p) and (p, ^2)- We repeat this process 
until no boundary vertices of degree two exist. This will ensure that the total 
number of boundary vertices over all subgraphs generated is 0{n/ ^/r) and 
this will also be a bound on the total number of super edges generated. 

The shortest path algorithm of [4] can easily be extended to deal with 
super edges in addition to the dense distance graphs: simply regard each 
super edge as a dense distance graph consisting of two vertices and one edge. 
Hence, the total running time for computing shortest paths in the first phase 
is 0{{n/y/r) log^n) = 0{n). 

Second phase We also face the problem with overlapping subgraphs when 
converting the shortest paths consisting of dense distance graph edges to 
shortest paths in G for the second phase of the algorithm. Klein's algo- 
rithm in] can report the underlying path in G corresponding to a dense 
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distance graph edge in time proportional to the length of the path (this 
requires that G has constant degree which we can assume without loss of 
generality). Hence, after the first phase we can obtain each of the shortest 
paths computed in time proportional to their total size. However, this size 
can be super-linear since the paths can share many vertices. We deal with 
this problem in the following. We will show that an implicit representation 
of the paths can be computed in 0{n) time. 

Implicit representation of paths Define pj^, . . . ,pj^ as above, i.e., the 
ordered sequence of vertices of P from which the coarse Reif algorithm has 
computed shortest paths Pi-^^, . . . ,Pi^. Let p^^ , . . . , p[^ be the endpoints of 
these paths. We start by obtaining the shortest path Qji in G from pi-^ to p^^ 
using Klein's algorithm on each dense distance graph edge of Pjj. This takes 
Oi\Qi,\) time. 

To find the shortest path Qi^ in G from pi^ to p[^, we similarly apply 
Klein's algorithm on Pi^. If we encounter no vertices already visited, we 
obtain the entire path and move on to Pi^. Otherwise, let t^i be the first 
already visited vertex and let Q^-^^ be the path found. We stop the algorithm 
when reaching Vi and instead start obtaining vertices of Qi^ backwards from 
p^2 until reaching an already visited vertex V2. Let Q^^ be the path found 
but ordered from t>2 to p[^. 

If V2 is on Qji , a simple property of shortest paths allows us to choose Qi^ 
as the concatenation of Q^^^, the subpath of Qi-^ from vi to V2, and Qy^, in 
that order. This gives us an implicit representation of Qi^ in O^Qi^ \ QiJ) 
time. 

In a dense distance graph, an edge [u, v) need not represent the same 
underlying shortest path as the edge (f , u) since shortest paths need not be 
unique. Hence, it may happen that V2 is on Q^^ and not on Qi^. If so, we 
can redefine Qi^ to be the subpath of Q^^ from pi^ to V2 followed by Q„2- 
Letting Qv2,vi be the subpath of Q^^ from V2 to fi, total time to find Qi^ 
is 0(|(5i2l + |Q»;2,t)il) = 0(|Qj2 \ Qn\ + \Qv2,vi\)- Since none of the vertices 
on Qv2,vii excluding V2, will be visited again, we can afford to spend time 

Repeating this process for the remaining shortest paths Pj3 , . . . , Pi^ gives 
an implicit representation of the corresponding shortest paths in G and it 
follows from the above analysis that running time is 0{n). In linear time it is 
then easy to obtain from this representation the desired subgraphs needed in 



the second phase of our algorithm and to ensure that they have total linear 
size. 

5 r-division 

In this section, we prove Theorem [H i.e., we show that an r-division can 
be found in 0{nlogr + (ra/i/r) logn) time. We use an approach similar to 
that of Frederickson [6]: contract G, find an r-division of this smaller graph, 
expand the graph back to G, and split some of the resulting pieces further to 
get the desired r-division of the whole graph. First, however, we shall give a 
simple 0(?7,log?T,) time algorithm. In Section ESI this algorithm will be used 
to find an r-division of the contracted graph. 

5.1 Weak r-division 

To obtain an r-division of G in 0{n log n) time, we again follow Frederickson 's 
approach. First we find a weak r-division which is a division of G into 0{n/r) 
pieces each of size 0{r) and with a constant number of holes, such that the 
total number of boundary vertices over all pieces is 0{n/^/r). In Section [5^ 
we will then split pieces further to get the desired r-division in O(nlogn) 
time. 

Consider the following recursive algorithm to find a weak r-division: re- 
gard G as a piece with no boundary vertices, split it recursively into two 
subpieces with the cycle separator theorem of Miller and recursive on them. 
The recursion stops when a piece has size at most r. As shown by Frederick- 
son [6], this gives a weak r-division. However, it does not ensure a constant 
bound on the number of holes in each piece which we need in our application. 

To deal with this, we use ideas of Fakcharoenphol and Rao [1] to keep the 
number of holes bounded by some constant h. The initial piece is the whole 
graph and thus contains no holes. Now, consider the general recursive step 
and let P = (Vp, Ep) be the current piece. Assume it has at most h holes. 
Apply the cycle separator theorem to P with all vertices assigned weight 
l/|Vp|. This splits P into two subpieces Pi = (Vp^, -EpJ and P[ = (Vp', Ep'), 
where |VpJ = a\Vp\+0{^/\V^\) and \Vp>\ = {l-a)\Vp\ + 0{./\V^\) for some 
I < ct ^ §• Assume w.l.o.g. that Pi belongs to the interior of the separator 
cycle. Then this subpiece has at most h holes. However, since the separator 
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cycle may have introduced a new hole in P{, this subpiece may have h + 1 
holes. 

Contract the holes of P[ into super vertices and apply the cycle separator 
theorem with vertex weights distributed evenly on super vertices. Expand 
them back to holes and let P2 and P3 be the resulting two pieces. As shown 
in [1], the number of holes in each of the two subpieces will be a constant 
factor smaller than /i + 1 so if we pick h sufficiently large, P2 and P3 each 
have at most h holes. Now, we recurse on Pi, P2, and P3 until all pieces 
contain at most r vertices. 

Lemma 2. The above procedure gives, for any parameter r G (0,-^,), a weak 
r-division of G where each piece has a constant number of holes. Running 
time is 0{n\og{n/r)). 

Proof. We have already argued that the number of holes in the pieces gener- 
ated is constant. We now show that the total number of boundary vertices 
over all pieces is 0{n/ ^). 

For any boundary vertex v in the weak r-division, let h{v) denote one less 
than the number of pieces containing v. Let B{n) be the sum of h{v) over all 
such V. There are nonnegative values ai, 02, and a^ with ai + 02 + 0^3 = 1 
such that in the above procedure, Pj contains at most a^n + c' y/n vertices, 
z = 1, 2, 3. Note that | < «i < |. We shall assume that a-z > 03. For n > r, 

B{n) < cy/n + B{ain + d^pn) + Bia^n + cy/n) + B{a^n + dy/n) 

for constants c, d > 0, and B{n) = for n < r. We will prove by induction 
on n > I that B{n) < d{n/y/r — ^y/n) for some constant d > {to be 
specified) . 

Clearly, this holds for ^ < n < r ioi any choice of c? > so assume that 
n > r and that the claim holds for smaller values. We have ai > | and since 
«2 > «3 and 02 + «3 > |, we have 02 > |- Hence, both aiu + dy/n and 
a2n + dsfn are at least ^ > | so the induction hypothesis can be applied 
to both of these values. We distinguish between two cases: a-^n + c'-y/n > | 
and a-^n + d ^fn < |. Assume first that a^n + d^/n > |. The induction 
hypothesis gives 



B{n,r) < cy/n ■ 



dn "iddy/n 



r \ r 



a\n + d\pn + \ a2n + c'y/n + \ a-iU + dy/n ) . 
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We will prove that the right-hand side is at most d{n/ y/r — \\/n), i.e., 
that 



Ud d d / d / d 



which will follow if we can show that 
. 3c 9c' 



, . ^ < y/al + ^/a^ + y/ai. (1) 

d Jr 



We may assume that r is at least some large constant. Picking d suffi- 
ciently large as well, we can make the left-hand side in ([1]) equal to 1 + e for an 
arbitrarily small constant e > 0. Since 02, «3 < 1 and since ai + 02 + ^3 = 1, 

a/oi + A/a2 + A/tt3 > (\/ai — ai) + ai + «2 + as = a/oi — cti + 1. 



Since | < ai < |, the right-hand side in ([T]) is larger than yf^l + l > 1- 

This proves the induction step for the case a^n + d ^/n > |. 

Now, assume that a^n + d^/n < |. Then B{a^n + d^/n) = and the 
induction hypothesis gives 

^ dn 2ddy/n d ( I ^ I ^\ 

5(n, r) < ci/n H — -= H ^ 77 V "i'^ + c v '^ + V (^2"^ + c'v'^ • 

Vr Vr 3 V^ V / 

The induction step will follow from the inequality 

3c 6c' , , . , 

1 + ^ + ^ < v^ + v^- (2) 

d Jr 



Since 

a^n < a^n + dy/n < - < — , 

we have 03 < | and hence ai + a2 > |. Since also | < «i < |, the right-hand 
side of ([2]) is greater than 

/ ^ N /2 2 8 ^ 

(Vtti -ai) + ai + a2>A/--- + ->l 

and ([2]) follows by picking d and r sufficiently large. 
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We have shown that the total number of boundary vertices over all pieces 
is 0{nj ^\ To show that the procedure generates a weak r-division, we also 
need to give an 0{n/r) bound on the number of pieces. Pieces of the form Pi 
or P2 each have size 0(r). Since the total number of vertices over all pieces 
is n + B{n) = 0{n), the number of such pieces is 0{n/r). The number of 
pieces of the form P3 cannot be larger than the number of the form Pi (or 
P2). Hence, the total number of pieces is 0{n/r). 

It follows that the procedure generates a weak r-division. Since Miller's 
cycle separator can be found in linear time and since the procedure recurses 
until pieces have size at most r, running time is 0{nlog{n/r)). D 

5.2 r-division in 0(n log n) time 

To obtain an r-division in 0{nlogn) time, we first find a weak r-division 
with Lemma El Each piece has 0{r) vertices and a constant number of holes 
but there may be more than order ^/r boundary vertices in the piece. 

We continue to follow Frederickson's approach while ensuring a constant 
number of holes in each piece. If there is a piece P containing more than c^/r 
boundary vertices for some constant c, apply the cycle separator theorem as 
in the weak r-division procedure to obtain subpieces Pi and Pi'. However, 
instead of distributing the vertex weights evenly on all vertices of P when 
applying the theorem, we now distribute weights evenly on boundary vertices 
only. We then infer subpieces P2 and P3 of P[ as before by distributing vertex 
weights evenly on super vertices defined by contracted holes of Pi'. This is 
repeated until each piece has at most c^/r boundary vertices. 

Lemma 3. The above procedure gives, for any parameter r G (0,n), an 
r-division of G in 0(n log n) time. 

Proof. The proof is more or less identical to that in [6] . We include it here 
for completeness. In the weak r-division, let tj be the number of pieces 
with exactly i boundary vertices. From the proof of Lemma [21 we have 
'^i^ti = S^gy (^('^) + 1)) where Vb is the set of boundary vertices over all 
pieces in the weak r-division. Hence, ^^gy iK'^) + 1) < '2B{n) = Oinj \fr). 
In the weak r-division, consider a piece P with i > c^ boundary vertices. 
When the above procedure splits P into subpieces Pi, P2, and P3, each 
of them contains at most a constant fraction of the boundary vertices of 
P. Hence, after di/{cy/r) splits of P for some constant rf, all subpieces 
will contain at most Cy/r boundary vertices. This will result in at most 
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1 + di/{cy/r) subpieces and at most c'a/t new boundary vertices per split for 
some constant c' . We may assume that c' < c. The total number of new 
boundary vertices introduced by the above procedure is thus 

Y,{c'^){di/{c^))U < dJ2^U = 0{n/^) 

i i 

and the number of new pieces is at most 

J2idt/icV^))t, = Oin/r). 

i 

Hence, the procedure generates an r-division. Since a weak r-division can 
be found in 0{n logn) time, an r-division can also be found within this time 
bound. n 

5.3 A faster algorithm 

We now show how to get the desired running time of 0(n \ogr + {n/ y/r) log n) 
in Theorem [T] We start by computing a spanning tree T of G (here, we as- 
sume that G is connected; we can always add infinite- weight edges to achieve 
this) and partitioning it into Q{n/ \/r) subtrees each of size 0(\/r); the sub- 
trees cover all vertices and are pairwise vertex-disjoint. This takes 0{n) time 
with the algorithm in [5]. 

Let G' be a plane multigraph obtained from G by contracting each subtree 
to a single vertex. This graph contains Oinj ^^ vertices. To obtain the same 
asymptotic bound on the number of edges, we will turn G' into a so called 
thin graph. In a plane multigraph, a higon is a face defined by two vertices 
and edges. A plane multigraph is thin if it contains no bigons. The following 
result from [1] shows that thin multigraphs are sparse. 

Lemma 4. A thin n-vertex multigraph contains 0{n) edges. 

Let G" be the thin multigraph obtained from G' by identifying the two 
edges of a bigon with one edge and repeating this process until no bigons 
exist. Graph G" is turned into a simple graph G'" by subdividing each edge 
{u, v) into two edges (u, w) and {w, v). Note that {u, v) corresponds to 0{y/r) 
edges in G' and in G; we subdivide each of them similarly such that the sum 
of weights of each edge pair equals the weight of the edge they subdivide. 
Now, G' is a simple graph and we colour black those vertices of G' that 
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correspond to contracted trees in G. All other vertices of G' are coloured 
white. 

By Lemma m G'" is a simple planar graph of size 0{n/^/r) so we can find 
an r-division of it in O {{n/^/r) log n) time with Lemma [31 This r-division 
consists of Oin/r^^"^) pieces each of size 0{r). We get an induced division 
of G' into pieces each consisting of 0{r) black vertices and 0{r^^'^) white 
vertices. Furthermore, each piece in this division has 0{^Jr) black boundary 
vertices and 0{r) white boundary vertices. 

Let Vi be the set of subtrees of G defined by the expanded black boundary 
vertices of pieces in the division of G' . Note that \Vi\ = 0{n/r) and each 
subgraph in Vi has size 0{y/r). 

For each piece P in the division of G", let P' be the piece in G defined by 
the union of edges in P and subtrees from expanded black interior vertices 
of P. Let V2 denote the set of these pieces P' . Note that IP2I = 0(n/r^/^) 
and each piece P' in V2 has size 0(r^/^). Furthermore, since there are 0{r) 
white boundary vertices in P' and each of the 0{y/r) black boundary vertices 
of the corresponding piece in G' contributes with at most 0{\/r) boundary 
vertices to P' when expanded, P' has 0{r) boundary vertices. 

The pieces in Vi U V2 cover all edges in G and each edge is contained in 
exactly one piece. We will transform these pieces into an r-division of G. 

First consider pieces P G "Pi. The number of boundary vertices of P is 
bounded by the size 0{y/r) of P. Since P is a tree, all its boundary vertices 
are on the external face. Hence, P has no holes and we include it as part of 
the r-division of G. 

Now, consider pieces P E V^- The piece P' in the division of G' cor- 
responding to P has a constant number of holes. We claim that the same 
holds for P. For consider some black boundary vertex v in P' and let Ty be 
the tree in G obtained by expanding v. In P, v gets expanded into 0{^/r) 
boundary vertices all belonging to T^. Since none of the edges of Ty belong 
to P by definition, these boundary vertices must all be on the same face of 
P; see Figure [2l Repeating this argument for all black boundary vertices of 
P', it follows that P has a constant number of holes. 

Since P has size 0{r^^'^) and 0{r) boundary vertices, we can find an r- 
division of it in 0{r^^'^ logr) time using Lemma [3] with a small modification: 
when finding a weak r-division of P, the boundary vertices and the holes of 
P will be regarded as boundary vertices and holes in the initial graph P. The 
result will still be a weak r-division of P since the total number of boundary 
vertices will be 0{\P\/^/¥ + r) = 0{r) = 0{\P\/^/¥). 
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Figure 2: (a): A piece P' in the r-division of G' with a black boundary 
vertex v on the external face, (b): After expanding f to a tree when forming 
a corresponding piece P G P2, new boundary vertices will also be on the 
external face. The same is true for holes. Only solid edges are part of the 
pieces. 

Total time to find r-divisions over all P G P2 is O(nlogr). Taking the 
union of the pieces obtained in all these r-divisions together with the pieces 
in Vi, we obtain the r-division of G in 0(nlogr + {n/^/r) logn) time. This 
proves Theorem [H 

6 Concluding Remarks 

We showed how to compute a min st-cut of a planar undirected n-vertex 
graph in 0(n log log n) time, improving on an earlier O(nlogn) bound. Can 
we get linear running time? Does a matching time bound also hold for planar 
directed graphs and for the maximum st-flow problem in planar (directed or 
undirected) graphs? 

The multiple source shortest path algorithm of Klein |9] has been general- 
ized to bounded genus graphs [3] . Fakcharoenphol and Rao |1] mention that 
their algorithm might apply to such graphs as well. It is therefore natural to 
conjecture that our time bound also holds for bounded genus graphs. 
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